Perturbative QCD study of ^(5) (f)p decays 

Jing-Wu W* Fu-Yi You^t 
Department of Physics, Xu Zhou Normal University, XuZhou 221116, China 

February 2, 2008 

o : 
o . 

(N ■ 

' Abstract 

We study -B(s) — > (f)p decays in a perturbative QCD approach based on kx factorization. In this 
, approach, we calculate factorizable and non-factorizable contributions, there are no annihilation 

' contributions due to quark content. We get the branching ratios and polarization fractions for 

T-H , 0p decays . Our predictions are consistent with the current experimental data. 
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tJ- ■ 1 Introduction 
(N ; 

■ Exclusive B meson decays, especially B VV modes, have aroused more and more interest both theo- 
\^ , retically and experimentally. Since the first observed charmless B VV mode, the B — > (j)K* decay [1], 
O I many B l/^decay channels have been studied in PQCD approach, such as i? ^ K* p, K*u![2], 
^ . B ^ K*K*[3], Bs p{co)K*[4], B p{uj)p{Lu)[5], and S° cl)(j)[6]. It offers an excellent place 
to study the CP violation and search for new physics hints [7]. Because the hadronization process 
Qh! is non-perturbative in nature, the essential problem in handling the decay processes is the separa- 
^ ' tion of different energy scales, i.e., the factorization assumption. Many approaches based on the 
factorization assumption have been developed, such as the naive factorization [8], the generalized 
• w^ ', factorization [9, 10], the QCD factorization[ll], and the perturbative QCD approach which is based 
^ on kx factorization[12, 13]. 

' Recently, B (pK* data reveal a large transverse polarization fraction, which has been considered 

as a puzzle, many theoretical efforts have been put to clarify it[14, 15, 16, 17, 18, 19, 20, 21]. This 
suggests that B — > VV modes must be more complicated than all the other modes and need to be 
studied deeply. Motivated by this, we study another B — > VV mode in the perturbative approach 
(PQCD) within the Standard Model. For B ^ (pp decay, only penguin operators contribute and we 
find that the branching ratio is at the order of 10~^. For B^ — > (fyp^ decay, current-current operators 
and penguin operators can contribute and we find that the branching ratio is at the order of 10~^. 
The longitudinal polarization predominates over transverse polarization and its fraction is found to 
go beyond 70%. Our predictions are consistent with the current experimental values. We hope that 
our study will help to resolve the above-mentioned puzzle a bit. 

The remaining part of this paper is organized as follows. In Sec. II, we calculate analytically the 
related Feynman diagrams and present the various decay amplitudes for the decay modes studied. In 
Sec. Ill, we give the numerical analysis for the branching ratios and polarization fraction of the related 
decay modes and compare them with the measured values. The summary and some discussion are 
included in the final section. 
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2 Theoretical framework and perturbative calculation 



In PQCD approach, the decay amplitude is expressed as the convolution of the mesons' light-cone 
wave functions, the hard scattering kernels and the Wilson coefficients, which stand for the soft, hard 
and harder dynamics, respectively. The formalism can be written as: 

~ j dxidx2dx3bidbib2db2b3db3Tr[C{t) 

6i)$<^(x2, b2)^p{x3, 63) 

H{xi,bi,t)St{xi)e-^^'^] (1) 

where Tr denotes the trace over Dirac and color indices. C(t) is Wilson coefficient of the four-quark 
operator which results from the radiative corrections at short distance. The wave function <^m absorbs 
non-perturbative dynamics of the process, which is process independent. The hard part H is rather 
process-independent and can be calculated in perturbative approach. The hi is the conjugate space 
coordinate of the transverse momentum, which represents the transverse interval of the meson, t is 
the largest energy scale in hard function H, while the jet function St{xi) comes from the resummation 
of the double logarithms In^Xj, called threshold resummation [22, 23], which becomes larger near the 
endpoint. The Sudakov form factor S{t) is from the resummation of double logarithms In^ Qb [24]. 

In this paper, we use the light-cone coordinates to describe the four-dimensional momentum as 
{p^ ,Pt)- We work in the frame with the B meson at rest, so the meson momentum can be written 
as 

P2 = ^{l-r%r%0T), 

Ps = ^irl,l-rl,OT) (2) 

in which r^, is defined by Vp = Mp/Ms and = M^/Mq. Pi, P2, P3 refer to B, cp, p respectively. To 
extract the helicity amplitudes, we parameterize the following polarization vectors. The longitudinal 
polarization must satisfy the orthogonality and normalization: £21 • -P2 = 0) ^3L ■ -P3 = 0) and = 
= — 1. Then we can give the manifest form as follows: 

As to the transverse polarization vectors, we can choose the simple form: 

£2T =^(0,0, It) 

esT =^(0,0, It) (4) 
The decay width for these channel is: 

r = i^ 1:m-M' (5) 

■fc* a=L,T 
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where |Pc| is the three-dimensional momentum of the final state meson, and |Pc| = "^(1 ~ ''p ~ f"^)- 
The subscript a denotes the hclicity states of the two vector mesons with L{T) standing for the 
longitudinal (transverse) component. As discussed in Ref.[l], the amplitude M'^ is decomposed into 

M" = mIMl + MlMNe*2{u = T) ■ el{a = T) 

+iMTSi,,pa4*4*P2Pi (6) 
We can define the longitudinal Hq, transverse H± helicity amplitudes 



Ho = mIMl, H± = mIMn T M^MpVr^ - IMt (7) 
where r = P2 • P^/{M^Mp). And we can deduce that they satisfy the relation 

M"^M" = \Ho\^ + \H+\^ + \H_\^ (8) 

a=L,T 

There is another set of definition of helicity amplitudes 

= ^M^Mp^2{r^ - 1)Mt (9) 



with the normalization factor ^ = y G|,Pc/(167rM^r). These helicity amplitudes satisfy the relation, 

|^o|' + |^|||' + |^±|' = l (10) 

where the notation Ao,^||,^_l denote the longitudinal, parallel, and perpendicular polarization am- 
plitudes. 

Our next task is to calculate the matrix elements M^jMn and Mt of the operators in the weak 
Hamiltonian with PQCD approach. We have to use the mesons' light-cone wave functions, they are 
universal for all decay channels. We employ the following wave functions as in other PQCD calculations 
[2, 3, 4]. 

1 

:/Pi+Mb)75$b(x,6) , 



V2A^ 

^[M, h{TWM)+ MT) m^pix) + ^^ie,,p,^,re^,{T)Pinl^;{x)] . 

where n+ = (1,0, Or) and n_ = (0, 1,0t) are dimensionless vectors on the light cone, x is the 
momentum fraction. 



3 



2.1 B ^ (j)p decays 

The effective Hamiltonian for the process^ — (j)p is given as [25] 

Heff = ^{Vu[c,{pmp) + C2{pmp)] 
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where Vu = VudYub ) = ^tdYtb ■> Ci{fJ,) are the Wilson coefficients, and the operators are 

Oi = {diUj)v-A{ujbi)v-A, O2 = {diUi)v-A{Ujbj)v-A, 

O3 = {dibi)v-A^{qjqj)v-A, O4 = {dibj)v-AY^{qjqi)v-A, 

O5 = {dibi)v-AJ2(^j'io)v+A, Oq = {dibj)v-A^{qjqi)v+A, 

1 Q 

3 - 3 - 

O7 = 2i^ibi)v-A^eq{qjqj)v+A, Os = -{dibj)v-A'^eq{qjqi)v+A, 
q q 

3 - 3 - 

O9 = ^{dibi)v-A'^eq{qjqj)v-A, Ow = ^{dibj)v-A'^eq{qjqi)v-A- 
q 9 



(11) 



(12) 



Here i and j stand for SU{3) color indices. The sum over q runs over the quark fields that are active 
at the scale p = 0(mi,), i.e., (q G {u,d,s,c,b}). From the effective Hamiltonian, we can see that the 
current-current operators have no contribution. For factorizable diagrams, all the penguin operators 
contribute, but for the non-factorizable diagrams only the operators 04,Oq,Os,Oio can contribute 
because of the color structure. The leading order diagrams for these decays are shown in Fig.l. We 
first calculate the usual factorization diagrams (a) and (b). 



Fu = SttCfM^ dxidxs bidbibadb^^sixiM) 

X{ [(1 + X3)^>p(x3) + - 2X3)($* (X3) + $^(X3)) 

xEe{t''e^)he{xi,X3,bi,b3) 
+2rp^i{x3)Ee{t^^^)he{x3, XI, 63, 61)} 



(13) 



Fjve = S-kCfM^ Jq dxidxs bidbib3db3^Bixi,bi 
xr^{ [$J(x3) + 2rp^>^(x3) + rpX3($^(x3) - $p(x3)) 

xEe{t''e^)he{xi,X3,bi,b3) 
+rpm{x3) + ^^Jx3)]Ee{t^e'')he{x3, X,, 63, ^l)} 



(14) 



Fxe = WttCfM'^ dxidx3 /q°° bidbib3db3^B{xi, bi 
xr^{[<^^{x3) + 2rp$;;(x3) - rpX3($^(x3) - '^p{x3)) 

xEe{t''P)he{xi,X3,bi,b3) 
+rp[^l{x3) + ^<i{x3)]ESf'^)he{x3, X^ 63, ^l)} 



(15) 
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where Ci? = | is a color factor. The function he, including the jet function S't(x3) (threshold resumma- 
tion for non-factorizable diagrams is weaker and negligible ), is the same as the kg in [1] . The factors 
E{t) contain the evolution from the W boson mass to the hard scales t in the Wilson coefficients a{t), 
and from t to the factorization scale 1/6 in the Sudakov factors S{t): 



Ee{t) = asit)aeit)SBit)Spit) 
The Wilson coefficients ae{t) in Eq.l6 are given by 

ae(t)-C3 + - + C5 + --- ------ 



(16) 



(17) 



For the non-factorizable diagrams (c) and (d) , all the three meson wave functions are involved and 
the amplitudes Mhc = M-Hei + M-HeG are written as 



MLe4 = IQ-KCpMl^/We dXidX2dX3 61^6162^62*5 (s^l , 
X{^^{X2) [-{X2 + X3)$p(x3) + rpX3($* (X3) + ^^(xs)) 

X £;e4 (4^^ ) (Xi , X2 , X3 , 61 , 62 ) 
+^^ix2) [(1 - X2)$p(.X3) + rp.X3($* (Xa) - «>^(X3)) 

X Ee4 {t^P ) (Xi , X2 , X3 , 61 , 62 ) 

MncA = IQirCpMly/Wc Q dxidx2dx3 hidhih2dh2'^B{xi, hi) 
xr^{[a;2(*;(x2) + ^4>''{x2))^1{x^) 
-2rp{x2 + X3)($^(x2)$;(a;3) + ^%{,X2)<l^^p{xM 

XEeiit^J^)h^J\xi,X2,X3,bi,b2) 
+ (1-X2)($;(X2) + $;(X2))<1>J(.X3) 

X Ee4 {t'^^^ )h'-^\xi, X2, X3,bi,b2)} 

MtsA = S2^^CFM^^/2N'c dxidx2dx3 6id6i62£i62$s(xi, 61) 

Xr^{[x2i<^l{x2) + $0"(X2))$J(X3) 
-2rpiX2 + X3)($;(X2)$«(X3) + $^(X2)$;;(X3))] 

X £'e4 (4^^ ) /i^^^ (Xi , X2 , 2:3 , 61 , 62 ) 
+ (1-X2)($^(X2) + $^(X2))#J(X3) 

XEei{t^^^)h^^^ (Xl, X2, X3, 61, 62)} 

•M.Le& = —'^Qi^C pM'^yJlNc /g^ dxidx2dx3 /g°° bidbib2db2^ b{xi, bi) 
x^^{x2){ \x2'i>p{x3) + rpa;3($* (xs) - ^pixs)) 



+ 



xEe6{tll^)hll\xi,X2,X3,bi,b2) 
-(1 - X2 + X3)$p(x3) + rpX3($* + %{X3)) 

X Ee6 (4^^ ) ^d^^ (a^l , X2 , X3 , 61 , 62 ) 



(18) 



(19) 



(20) 



(21) 
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MncG = -IQirCpM^^JlNc dxidx2dx^ bidbib2db2^B{xi, bi) 
xr^{x2{%{x2) - n^{x2))'^^{x3)Eee{t'-J^)h^J\xi,X2,X3,bi,b2) 

+ [(1-X2)(«J>;(X2)-<I>;(.T2))$^(X3) 
-2rp(l -X2+ X3)(«J>^(X2)$^(X3) - n(^2)'^>M 

xE^Q{tf)hf (Xi, X2, X3, hi, 62)} (22) 

MTe& = -327rCFM2 V2]Ve /o dxidx2dx3 6id6i62d62«'i?(xi, 61) 
xr4x2($;(x2) - $(/.'^(x2))$^(a;3)^e6(4'V?^(a^i'^2,X3,6i,6 2) 

+ [(1-X2)(<&;(X2)-<I>;(X2))<I>^(X3) 

-2rp(i - X2 + X3)(«5;(x2)$^(x3) - $;(x2)^>;;(x3))] 

X£;e6(4'^)/lf (Xi, X2, X3, 61, 62)} (23) 

The evolution factors are given by 

Eei{t) = asit)ai{t)Sit)\ 

63=61 

(24) 

with the Sudakov factor S = SsS^Sp. The Wilson coefficients a appearing in the above formulas are 

"4—3 6 ' 

«6 = f - f (25) 
The amplitudes forS+ — are written as 

Mh = UV^FHe + V^Mue (26) 

where Fhc denotes the factorizable contributions and Mue the non-factorizable contributions. For 
the other two decay channels of 5 — > (j)p, the amplitudes are the following: for B~ (pp~ , 

Mh = UVtFHe + VtMHe (27) 

and for 5° ^ (j)p°, 

-V2M% = V*UFHe + Vt*MHe (28) 

2.2 5° ^ decay 

For the process Bg (pp^ , it is a 6 — > s transition and we use the effective Hamiltonian [25] 

^ 10 
-Heff = [VubV:,{CiO^ + C2O2") - Vt,VZ (^ a Oj"))] . (29) 

i=3 

We specify below the operators in T^g^^ for 6 ^ s: 

O" = {siUj)v-A{ujbi)v-A, O2 = {siUi)v-Aiujbj)v-A, 
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Os = {sibi)v-A^{qjqj)v-A, O4 = {sibj)v-A^{qjqi)v-A, 

q q 

O5 = {sibi)v-A^{qjqj)v+A, Oq = {sibj)v-A^{qjqi)v+A, 

3 3 
O7 = 2(^ibi)v-AY^eq{qjqj)v+A, Og = -{sibj)v-A'^eq{qjqi)v+A, 
q 1 

3 3 
O9 = 2(^ibi)v-A^eq{qjqj)v-A, Oio = ^isibj)v-A'^eq{qjqi)v-A- (30) 

q q 

Prom the effective Hamiltonian we can see that the current-current operators and penguin operators 
can contribute. The leading order diagrams for the decay are shown in Fig.2, the amphtude for 
Bg — p^<p mode is 

-V2M'h = V:f,F'H, + V:M'He - VCfpFj-', - v;Mfr'e (31) 

where M.^'^ = M.^^^ + M.^^^. Here Vu = Kib^s; = ^fe^ts and amphtude for the corresponding 
CP conjugate model is written as 

-^M'h = VufpF'H, + VuM'He - VtfpFZ - VtM^il, (32) 
Next we calculate the hard part in PQCD approach. Por the factorizable diagram (e) and (f), we have 



X{[(1 + X2)^^{X2) + r;(l - 2X2){%{X2) + %{X2)) 

xE'Ati^^)he{xi,X2,bi,b2) 
+2r'^^l{x2)E',{tP)he{x2,xi,b2M)} 

F'j,^ = SttCfMI^ dxidX2 /o°° bidbib2db2^BAxi, 6i) 
xr'p{ [^1{X2) + 2r;$-(x2) + r^X2{%{x2) - %{x2)) 

xE'At^e^)he{xi,X2,bi,b2) 
+^;[^;(^2) + %{x2)]EUt''^^)he{x2,Xl,b2,bl)} 

F^^ = mTTCpMl dXidx-i bidbib2db2^BAxi, bl) 

xr',{[<^>l{x2) + 2r;$;(x2) - r'^X2i%{x2) - %{x2)) 

xE'At''P)he{xi,X2MM) 

+r'A%{x2) + %{x2)]Ee{t'i^)K{x2,xiM, bl)] 



(33) 



(34) 



(35) 



with r'^,r'^ equal to r^,rp except for Mb replaced by Mb^- The expression for F^'^ are the same as 
^He ^^^^ Wilson coefficients a replaced by . As before, the factor E'{t) are given by 



Efy{t) = asaf\t)SB{t)S^{t) 
where the Wilson coefficients are the following 



(36) 
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a = Ci + 2/3C2, 
= 3/2C7 + 1/2C8 + 3/2CC9 + I/2C10. 
For non-factorizable diagram (g) and (h), we find that 



(37) 



M'l,^ = WttCfM^^V^cIo dxidX2dX3 /o°° hdbib3db3 

^s.(a;i,6i) X {$p(a;3)[-(x2 +X3)$0(x2) + r^X2($^(a;2) 
+%{X2))] X Ef{t'^j^)h'i\x,,X2,xsM,b3) 

+^p{x3) [(1 - X3)^4x2)+r'^X2{%{x2) - $^(^2)) 
xE(f(if )/lf (xi,X2,X3,6l,63)} 

•^Le6 = -16^CfM2 V27v; dxidx2dx3 /o°° hidhih^dhz 
^BsixiM) X $p(.X3){[.T3$</,(.X2) +r^a;2($^(x2) - 

$i(x2))] X (iW)4'^(^l,^2,X3,6l,63) 



+ 



-(1 - X3 + X2)$^(X2) + r;^a;2($*/ar2) + $^(^2)) 

XE^^\tf)hf{x^,X2,X3MM) 



(38) 



(39) 



•^?^e4 = leTrC^Ml^ V2iVe/o dxidx2dxs bidbihdh 
^Bs{xi,bi) X r;,{[x3($^(a:3) + $;5(x3))$J(.T2) 

-2r;^(x2 + X3)(<5^(.T3)$;(X2) + ^1{X3)^1{X2))] 

X£;(f (4^))4^^(xi,X2,X3,6i,63) 
+ (1-X3)($;;(X3) + $«(X3))$J(X2) 

X£;(f (4'))4')(X1, X2, X3, 61, 63)} (40) 



•^Wee = -leTrCpMl ^27v;/o^ dxidx2dx3 j;^ bidbi 

b3dbs<^BAxi,h) X r;{[x3 ($;;(X3) - $^(X3)) 
$^(X2)] X £;(f (4^VW(xi,X2,X3,6l,63) 

+[(1 - X3) ($^(X3) - $«(X3)) $J(X2) - 2r; 

(1 - X3 + X2) ($^(X3)$^(X2) - $«(X3)«J>;(X2))] 

xi?(f(tf )42)(xi,X2,X3,6l,63)} (41) 



M'Te4 = 327rCFM|V2iVc Jo rfxicZx2dx3 /o°° 61^6163 

db^^BMM) X r;{[x3($^(x3) + $^(X3))$^(X2) 
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XEi''^'{t^J^)h'j\xi,X2,Xs,h,h) 

+ (l-X3)($;;(x3) + $^(ar3))$^(x2) 
X4'>(if (xi,X2,X3,6i,63)} (42) 



X^gg = -167rCFM|x/2iVc/o dx^dx^dx^ hdbi 

bsdh^BM, h) X r;{x3[$^(X3) - $»(X3)]$J(X2) 
X£;^«)'(4'^)4')(X1,X2,X3,61,63) 

+[(1 - X3) {^;{xs) - $«(X3)) $J(X2) - 2r; 

(1 - X3 + X2) ($;^(X3)$;(X2) - $^(X3)^>;(X2))] 

x4'^)'(if )/if (xi, X2, X3, 6i, 63)} (43) 

The evolution factors are given by 

E^^'{t)=as{t)af{t)Sm,=,, (44) 

with the Sudakov factor S = SBsS(j>Sp and the Wilson coefficients are given by 

a'l = C2/N,, 
< = 3/2C10/N,, 

a'e = 3/2Cs/Nc. (45) 

2.3 Numerical analysis 

The parameters used in our calculations are: the Fermi couphng constant Gp = 1.16639 x 10~^GeV~'^, 
the meson masses Mb = 5.28GeV, Mbs = 5.37GeF, Mp = 0.77GeV, = IMGeV, the decay 
constant fp = 0.205GeV, fj = 0.155Gey, = 0.237^6^", /J = Q.22QGeV , the central value of the 
CKM matrix elements 7 = 60°, \Vtd\ = 0.0075, \Vtb\ = 0.9992, = 0.0047, \Vus\ = 0.2196 and the 
meson lifetimer^ = l.Qbps,TBs — 1.461ps[26]. 

Using the above parameters, we get the branching ratios and helicity amplitudes of <pp 
decays (the helicity amplitudes are in Table 1) 

Br{B^ (j)p'^) = 4.1 X 10-^ 
Br(B° ^ (/)/9°) = 1.9 X 10-^ 
Br{B'^ <j)p°) = 3.09 X 10"^, 

Br{Bs° (t)p°) = 3.60 X 10-'^ (46) 

Compared with the averaged results of QCDF[27] 

BR{B- ^ p-(t)) = 5.5 X 10-^ 
BR{B^ (/)p°) = 2.5 X 10-9. (47) 

or those from naive factorization [28] 

Bi?(5° ^ p^(j)) = 2.92 X 10-^ (48) 
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our predictions for B ^ (f)p are consistent with those of QCDF, the predictions for S° p^cf) are 
consistent with the result of naive factorization, because the nonfactorization effects in Bg —>■ p^(f) are 
httle. 

Presently, only the experimental upper limits are available at the 90% confidence level [26] , 

Br{B+ (t)p+) < 1.6 X 10~^ 
Br{B^ 4>p^) < 1.3 X 10-^ 

5^(5° ^ (l)p°) < 6.17 X lO-''. (49) 

Obviously, our results are consistent with the data. Our predictions will be tested by the oncoming 
measurements. 

For B ^ (pp decays, only penguin operators contribute, so there is no direct CP violation in the 
B ^ (pp decays. 

For Bg (j)p^ decays, the CP asymmetry is time dependent 

Acp{t) - ^cp cos(Amt) + a,+£/ sin(Amt), (50) 

where Am is the mass difference of the two mass eigenstates of neutral Bg mesons. 
The direct CP violation parameter is defined 

_ \M\^-\M\^ 
"^^^ - \MV + \Mr ^ ^ 

The direct CP violation parameter we can get in B^ (pp^ is 

Afp{Bl ^ (t^p^^) =-8.0%. (52) 

Prom Table 1, we can find the longitudinal fraction in B^ (pp^ decays go beyond 70%, but the 
longitudinal fractions in B ^ <pp decays are very small, which is similar to B^ — p^p'^ decay mode [29]. 

In B ^ (pp decays, Oi 2m "Heff don't contribute via factorizable diagrams, the penguin operators 
contributing via factorizable diagrams are color suppressed, so that the nonfactorizable effects are 
the same order as the factorizable ones, which cause the B —>■ (pp decays not to be factorization 
dominated, besides the above reason, the nonfactorizable longitudinal amplitude is opposite in sign 
to that of the factorizable part, therefore the longitudinal amplitude gets a large cancellation between 
the factorizable effects and the non factorizable parts such that \Aq\'^ is reduced much. 

According to Ref.[30], we can get the B^g-j ^ p,(p vector meson transition form factors , 

B^ p, F(0) = 0.303, Ao(0) = 0.308, 
Ai(0) = 0.233,^2(0) = 0.208, 
Bg 0, 1/(0) = 0.430, Ao(0) = 0.363, 
^i(O) = 0.304,^2(0) = 0.276, 
which are consistent with the results with light cone sum rules [31]. 



3 Summary 

In this paper, we calculate the branching ratios and polarization fractions of B ^ (/)p and Bg 
(pp^ decays in perturbative QCD approach, the predicted branching ratios are compared with the 
experimental data and results obtained with other approaches. CP parameters in B^ (pp^ are given 
in our paper, we compared with the experimental values, our results are consistent with the current 
experimental data. 

We thank Dong-Sheng Du and Mao-Zhi Yang for helpful discussions and communications. 
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Channel 






1 -Aperp 1 ' 




0.14 


0.41 


0.45 




0.14 


0.41 


0.45 




0.78 


0.12 


0.10 




0.83 


0.09 


0.08 



Table 1: Branching ratios and helicity amplitudes 

A wave functions 

The (f) and p distribution amplitudes up to twist 3 are given by 





3/</. ... 




fT 


nix) 


fT 

- / [3(1 
4V2iVc 


^lix) 


3/J , 
= — --^ — xil 
2^J2Nc 


ni^) 


- ^ i \1 


ni^) 


4V2iV/ 



3(1 - 2x)(4.5 - 11.2a; + 11.2a;^) + 1.38 In 



1 + O.2C42 (1 - 2x) 



with the Gegenbauer polynomials, 



clio = lis5e-soe+^), 
clio = li5e-i). 



(53) 



(54) 



For p meson , its Lorentz structures are similar to (p meson and the distribution amplitudes are 

given by 



%ix) = 



Kix) 



nix) 



jfp_ 



2^/2N, 



xil -x)\l + 0.18C2^^(1 - 2x) 



{3(1 - 2xf + 0.3(1 - 2x)^[5(l - 2xf - 3] + 0.21[3 - 30(1 - 2xf + 35(1 - 2x)% , 



^ (1 -2x)[l + 0.76(10x2 - lOx + 1)] , 



3/J 



nix) = ^^xil -x)[l + Q.2Cl'\l - 2x) 
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=(1 - 2x) \l + 0.93(10x2 - lOx + 1) 



For the amplitudes of B andS^, meson, we employ the following distribution amplitudes: 



which satisfies the normalization 



Jo 



f. 



B 



(55) 



(56) 



We choose Nb = 91.784GeV, cob = 0.4GeV. Things for Bg are similar if we ignore SU(3) symmetry 
breaking effect. As discussed in Ref. [4], we choose cob, = O.SOGeV. 
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' Abstract 

We study -B(s) — > (f)p decays in a perturbative QCD approach based on kx factorization. In this 
, approach, we calculate factorizable and non-factorizable contributions, there are no annihilation 

' contributions due to quark content. We get the branching ratios and polarization fractions for 

B{s) 4'P decays . Our predictions are consistent with the current experimental data. 

(N 
> 

tJ- ■ 1 Introduction 
(N 

■ Exclusive B meson decays, especially B VV modes, have aroused more and more interest both 
\^ . theoretically and experimentally. Since the first observed charmless B — > VV mode, the B 

[ 4>K* decay [?], many B VVdecay channels have been studied in PQCD approach, such as i? — > 
^ ■ K*p,K*u[?], B K*K*[7], Bs p{uj)K*[7], B p{uj)p{uj)[7], and B^ It offers an excel- 

^p^. lent place to study the CP violation and search for new physics hints [?]. Because the hadronization 
Qh! process is non-perturbative in nature, the essential problem in handling the decay processes is the 
^ ' separation of different energy scales, i.e., the factorization assumption. Many approaches based on 
the factorization assumption have been developed, such as the naive factorization [?], the generalized 
factorization [?, ?], the QCD factorization [?], and the perturbative QCD approach which is based on 
^ ' kT factorization [?, ?]. 

' Recently, B (pK* data reveal a large transverse polarization fraction, which has been considered 

as a puzzle, many theoretical efforts have been put to clarify it[?, ?, ?, ?, ?, ?, ?, ?]. This suggests that 
B — > VV modes must be more complicated than all the other modes and need to be studied deeply. 
Motivated by this, we study another B — > VV mode in the perturbative approach (PQCD) within 
the Standard Model. For B ^ (pp decay, only penguin operators contribute and we find that the 
branching ratio is at the order of 10~^. For B^ — > (pp^ decay, current-current operators and penguin 
operators can contribute and we find that the branching ratio is at the order of 10"''. The longitudinal 
polarization predominates over transverse polarization and its fraction is found to go beyond 70%. 
Our predictions are consistent with the current experimental values. We hope that our study will help 
to resolve the above-mentioned puzzle a bit. 

The remaining part of this paper is organized as follows. In Sec. II, we calculate analytically the 
related Feynman diagrams and present the various decay amplitudes for the decay modes studied. In 
Sec. Ill, we give the numerical analysis for the branching ratios and polarization fraction of the related 
decay modes and compare them with the measured values. The summary and some discussion are 
included in the final section. 
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2 Theoretical framework and perturbative calculation 



In PQCD approach, the decay amplitude is expressed as the convolution of the mesons' light-cone 
wave functions, the hard scattering kernels and the Wilson coefficients, which stand for the soft, hard 
and harder dynamics, respectively. The formalism can be written as: 

~ j dxidx2dx3bidbib2db2b3db3Tr[C{t) 

6i)$<^(x2, b2)^p{x3, 63) 

H{xi,bi,t)St{xi)e-^^'^] (1) 

where Tr denotes the trace over Dirac and color indices. C(t) is Wilson coefficient of the four-quark 
operator which results from the radiative corrections at short distance. The wave function <^m absorbs 
non-perturbative dynamics of the process, which is process independent. The hard part H is rather 
process-independent and can be calculated in perturbative approach. The hi is the conjugate space 
coordinate of the transverse momentum, which represents the transverse interval of the meson, t is 
the largest energy scale in hard function H, while the jet function St{xi) comes from the resummation 
of the double logarithms In^Xj, called threshold resummation[?, ?], which becomes larger near the 
endpoint. The Sudakov form factor S{t) is from the resummation of double logarithms In^ Qb [?]. 

In this paper, we use the light-cone coordinates to describe the four-dimensional momentum as 
{p^ ,Pt)- We work in the frame with the B meson at rest, so the meson momentum can be written 
as 

P2 = ^{l-r%r%0T), 

Ps = ^irl,l-rl,OT) (2) 

in which r^, is defined by Vp = Mp/Ms and = M^/Mq. Pi, P2, P3 refer to B, cp, p respectively. To 
extract the helicity amplitudes, we parameterize the following polarization vectors. The longitudinal 
polarization must satisfy the orthogonality and normalization: £21 • -P2 = 0) ^3L ■ -P3 = 0) and = 
= — 1. Then we can give the manifest form as follows: 

As to the transverse polarization vectors, we can choose the simple form: 

£2T =^(0,0, It) 

esT =^(0,0, It) (4) 
The decay width for these channel is: 

r = i^ 1:m-M' (5) 

■fc* a=L,T 
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where |Pc| is the three-dimensional momentum of the final state meson, and |Pc| = "^(1 ~ ''p ~ f"^)- 
The subscript a denotes the hclicity states of the two vector mesons with L{T) standing for the 
longitudinal (transverse) component. As discussed in Ref.[?], the amplitude is decomposed into 

M" = mIMl + MlMNe*2{u = T) ■ el{a = T) 

+iMTSi,,pa4*4*P2Pi (6) 
We can define the longitudinal Hq, transverse H± helicity amplitudes 



Ho = mIMl, H± = mIMn T M^MpVr^ - IMt (7) 
where r = P2 • P^/{M^Mp). And we can deduce that they satisfy the relation 

M"^M" = \Ho\^ + \H+\^ + \H_\^ (8) 

a=L,T 

There is another set of definition of helicity amplitudes 

= ^M^Mp^2{r^ - 1)Mt (9) 



with the normalization factor ^ = y G|,Pc/(167rM^r). These helicity amplitudes satisfy the relation, 

|^o|' + |^|||' + |^±|' = l (10) 

where the notation Ao,^||,^_l denote the longitudinal, parallel, and perpendicular polarization am- 
plitudes. 

Our next task is to calculate the matrix elements M^jMn and Mt of the operators in the weak 
Hamiltonian with PQCD approach. We have to use the mesons' light-cone wave functions, they are 
universal for all decay channels. We employ the following wave functions as in other PQCD calculations 
[?, ?, ?]• 

1 

:/Pi+Mb)75$b(x,6) , 



V2A^ 

where n+ = (1,0, Or) and n_ = (0, 1,0t) are dimensionless vectors on the light cone, x is the 
momentum fraction. 
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2.1 B ^ (j)p decays 

The effective Hamiltonian for the process^ — (j)p is given as [?] 

Heff = ^{Vu[c,{pmp) + C2{pmp)] 
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where Vu = VudYub ) = ^tdYtb ■> Ci{fJ,) are the Wilson coefficients, and the operators are 

Oi = {diUj)v-A{ujbi)v-A, O2 = {diUi)v-A{Ujbj)v-A, 

O3 = {dibi)v-A^{qjqj)v-A, O4 = {dibj)v-AY^{qjqi)v-A, 

O5 = {dibi)v-AJ2(^j'io)v+A, Oq = {dibj)v-A^{qjqi)v+A, 

1 Q 

3 - 3 - 

O7 = 2i^ibi)v-A^eq{qjqj)v+A, Os = -{dibj)v-A'^eq{qjqi)v+A, 
q q 

3 - 3 - 

O9 = ^{dibi)v-A'^eq{qjqj)v-A, Ow = ^{dibj)v-A'^eq{qjqi)v-A- 
q 9 



(11) 



(12) 



Here i and j stand for SU{3) color indices. The sum over q runs over the quark fields that are active 
at the scale p = 0(mi,), i.e., (q G {u,d,s,c,b}). From the effective Hamiltonian, we can see that the 
current-current operators have no contribution. For factorizable diagrams, all the penguin operators 
contribute, but for the non-factorizable diagrams only the operators 04,Oq,Os,Oio can contribute 
because of the color structure. The leading order diagrams for these decays are shown in Fig.l. We 
first calculate the usual factorization diagrams (a) and (b). 



Fu = SttCfM^ dxidxs bidbibadb^^sixiM) 

X{ [(1 + X3)^>p(x3) + - 2X3)($* (X3) + $^(X3)) 

xEe{t''e^)he{xi,X3,bi,b3) 
+2rp^i{x3)Ee{t^^^)he{x3, XI, 63, 61)} 



(13) 



Fjve = S-kCfM^ Jq dxidxs bidbib3db3^Bixi,bi 
xr^{ [$J(x3) + 2rp^>^(x3) + rpX3($^(x3) - $p(x3)) 

xEe{t''e^)he{xi,X3,bi,b3) 
+rpm{x3) + ^^Jx3)]Ee{t^e'')he{x3, X,, 63, ^l)} 



(14) 



Fxe = WttCfM'^ dxidx3 /q°° bidbib3db3^B{xi, bi 
xr^{[<^^{x3) + 2rp$;;(x3) - rpX3($^(x3) - '^p{x3)) 

xEe{t''P)he{xi,X3,bi,b3) 
+rp[^l{x3) + ^<i{x3)]ESf'^)he{x3, X^ 63, ^l)} 



(15) 
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where Cp = | is a color factor. The function he-, including the jet function S't(x3) (threshold resumma- 
tion for non-factorizable diagrams is weaker and negligible ), is the same as the in [?] . The factors 
E{t) contain the evolution from the W boson mass to the hard scales t in the Wilson coefficients a{t), 
and from t to the factorization scale 1/6 in the Sudakov factors S{t): 



Ee{t) = asit)aeit)SBit)Spit) 
The Wilson coefficients ae{t) in Eq.?? are given by 

ae(t)-C3 + - + C5 + --- ------ 



(16) 



(17) 



For the non-factorizable diagrams (c) and (d) , all the three meson wave functions are involved and 
the amplitudes Mhc = M-Hei + M-HeG are written as 



MLe4 = IQ-KCpMl^/We dXidX2dX3 61^6162^62*5 (s^l , 
X{^^{X2) [-{X2 + X3)$p(x3) + rpX3($* (X3) + ^^(xs)) 

X £;e4 (4^^ ) (Xi , X2 , X3 , 61 , 62 ) 
+^^ix2) [(1 - X2)^pix3) + rp.X3($* (Xa) - «>^(X3)) 

X Ee4 {t^P ) (Xi , X2 , X3 , 61 , 62 ) 

MncA = IQirCpMly/Wc Q dxidx2dx3 hidhih2dh2'^B{xi, hi) 
xr^{[a;2(*;(x2) + ^4>''{x2))^1{x^) 
-2rp{x2 + X3)($^(x2)$;(a;3) + ^%{,X2)<l^^p{xM 

XEeiit^J^)h^J\xi,X2,X3,bi,b2) 
+ (1-X2)($;(X2) + $;(X2))<1>J(.X3) 

X Ee4 {t^^'')h''p {xi,X2,X3,biM)} 

MtsA = S2^^CFM^^/2N'c dxidx2dx3 bidhib2dh2^B{xiM) 

Xr^{[x2{^l{x2) + $0"(X2))$J(X3) 
-2rp{X2 + X3)($;(X2)$«(X3) + ^%{x2Wp{x3))] 

XEe4,it'j^)h^j\xi, X2, X3,bi,b2) 
+ (1-X2)($^(X2) + $^(X2))*J(X3) 

XEei{t^^^)h^^^ (Xl, X2, X3, 61, 62)} 

•M.Le& = —'^Qi^C pM'^yJlNc /g^ dxidx2dx3 /g°° bidbib2db2^ b{xi, bi) 
X^^{x2){ \x2'i>p{x3) + rpa;3($* (xs) - $^(^3)) 



+ 



X ^e6 (ii^"* ) /^i^^ (x 1 , a;2 , X3 , 61 , 62 ) 

-(1 -X2 + X3)$p(x3) + rpX3($* + %{X3)) 
XEe(i{t^^'')h^^\xi,X2,X3,hiM) 



(18) 



(19) 



(20) 



(21) 
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MncG = -IQirCpM^^JlNc dxidx2dx^ bidbib2db2^B{xi, bi) 
xr^{x2{%{x2) - n^{x2))'^^{x3)Eee{t'-J^)h^J\xi,X2,X3,bi,b2) 

+ [(1-X2)(«J>;(X2)-<I>;(.T2))$^(X3) 
-2rp(l -X2+ X3)(«J>^(X2)$^(X3) - n(^2)'^>M 

xE^Q{tf)hf (Xi, X2, X3, hi, 62)} (22) 

MTe& = -327rCFM2 V2]Ve /o dxidx2dx3 6id6i62d62«'i?(xi, 61) 
xr4x2($;(x2) - $(/.'^(x2))$^(a;3)^e6(4'V?^(a^i'^2,X3,6i,6 2) 

+ [(1-X2)(<&;(X2)-<I>;(X2))<I>^(X3) 

-2rp(i - X2 + X3)(«5;(x2)$^(x3) - $;(x2)^>;;(x3))] 

X£;e6(4'^)/lf (Xi, X2, X3, 61, 62)} (23) 

The evolution factors are given by 

Eei{t) = asit)ai{t)Sit)\ 

63=61 

(24) 

with the Sudakov factor S = SsS^Sp. The Wilson coefficients a appearing in the above formulas are 

"4—3 6 ' 

«6 = f - f (25) 
The amplitudes forS+ — are written as 

Mh = UV^FHe + V^Mue (26) 

where Fhc denotes the factorizable contributions and Mue the non-factorizable contributions. For 
the other two decay channels of 5 — > (j)p, the amplitudes are the following: for B~ (pp~ , 

Mh = UVtFHe + VtMHe (27) 

and for 5° ^ (j)p°, 

-V2M% = V*UFHe + Vt*MHe (28) 

2.2 5° ^ decay 

For the process Bg (pp^ , it is a 6 — > s transition and we use the effective Hamiltonian [?] 

^ 10 
-Heff = [VubV:,{CiO^ + C2O2") - Vt,VZ (^ a Oj"))] . (29) 

i=3 

We specify below the operators in T^g^^ for 6 ^ s: 

O" = {siUj)v-A{ujbi)v-A, O2 = {siUi)v-Aiujbj)v-A, 
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Os = {sibi)v-A^{qjqj)v-A, O4 = {sibj)v-A^{qjqi)v-A, 

q q 

O5 = {sibi)v-A^{qjqj)v+A, Oq = {sibj)v-A^{qjqi)v+A, 

3 3 
O7 = 2(^ibi)v-AY^eq{qjqj)v+A, Og = -{sibj)v-A'^eq{qjqi)v+A, 
q 1 

3 3 
O9 = 2(^ibi)v-A^eq{qjqj)v-A, Oio = ^isibj)v-A'^eq{qjqi)v-A- (30) 

q q 

Prom the effective Hamiltonian we can see that the current-current operators and penguin operators 
can contribute. The leading order diagrams for the decay are shown in Fig.2, the amphtude for 
Bg — p^<p mode is 

-V2M'h = V:f,F'H, + V:M'He - VCfpFj-', - v;Mfr'e (31) 

where M.^'^ = M.^^^ + M.^^^. Here Vu = Kib^s; = ^fe^ts and amphtude for the corresponding 
CP conjugate model is written as 

-^M'h = VufpF'H, + VuM'He - VtfpFZ - VtM^il, (32) 
Next we calculate the hard part in PQCD approach. Por the factorizable diagram (e) and (f), we have 



X{[(1 + X2)^^{X2) + r;(l - 2X2){%{X2) + %{X2)) 

xE'Ati^^)he{xi,X2,bi,b2) 
+2r'^^l{x2)E',{tP)he{x2,xi,b2M)} 

F'j,^ = SttCfMI^ dxidX2 /o°° bidbib2db2^BAxi, 6i) 
xr'p{ [^1{X2) + 2r;$-(x2) + r^X2{%{x2) - %{x2)) 

xE'At^e^)he{xi,X2,bi,b2) 
+^;[^;(^2) + %{x2)]EUt''^^)he{x2,Xl,b2,bl)} 

F^^ = mTTCpMl dXidx-i bidbib2db2^BAxi, bl) 

xr',{[<^>l{x2) + 2r;$;(x2) - r'^X2i%{x2) - %{x2)) 

xE'At''P)he{xi,X2MM) 

+r'A%{x2) + %{x2)]Ee{t'i^)K{x2,xiM, bl)] 



(33) 



(34) 



(35) 



with r'^,r'^ equal to r^,rp except for Mb replaced by Mb^- The expression for F^'^ are the same as 
^He ^^^^ Wilson coefficients a replaced by . As before, the factor E'{t) are given by 



Efy{t) = asaf\t)SB{t)S^{t) 
where the Wilson coefficients are the following 



(36) 
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a = Ci + 2/3C2, 
= 3/2C7 + 1/2C8 + 3/2CC9 + I/2C10. 
For non-factorizable diagram (g) and (h), we find that 



(37) 



M'l,^ = WttCfM^^V^cIo dxidX2dX3 /o°° hdbib3db3 

^s.(a;i,6i) X {$p(a;3)[-(x2 +X3)$0(x2) + r^X2($^(a;2) 
+%{X2))] X Ef{t'^j^)h'i\x,,X2,xsM,b3) 

+^p{x3) [(1 - X3)^4x2)+r'^X2{%{x2) - $^(^2)) 
xE(f(if )/lf (xi,X2,X3,6l,63)} 

•^Le6 = -16^CfM2 V27v; dxidx2dx3 /o°° hidhih^dhz 
^BsixiM) X $p(.X3){[.T3$</,(.X2) +r^a;2($^(x2) - 

$i(x2))] X (iW)4'^(^l,^2,X3,6l,63) 



+ 



-(1 - X3 + X2)$^(X2) + r;^a;2($*/ar2) + $^(^2)) 

XE^^\tf)hf{x^,X2,X3MM) 



(38) 



(39) 



•^?^e4 = leTrC^Ml^ V2iVe/o dxidx2dxs bidbihdh 
^Bs{xi,bi) X r;,{[x3($^(a:3) + $;5(x3))$J(.T2) 

-2r;^(x2 + X3)(<5^(.T3)$;(X2) + ^1{X3)^1{X2))] 

X£;(f (4^))4^^(xi,X2,X3,6i,63) 
+ (1-X3)($;;(X3) + $«(X3))$J(X2) 

X£;(f (4'))4')(X1, X2, X3, 61, 63)} (40) 



•^Wee = -leTrCpMl ^27v;/o^ dxidx2dx3 j;^ bidbi 

b3dbs<^BAxi,h) X r;{[x3 ($;;(X3) - $^(X3)) 
$^(X2)] X £;(f (4^VW(xi,X2,X3,6l,63) 

+[(1 - X3) ($^(X3) - $«(X3)) $J(X2) - 2r; 

(1 - X3 + X2) ($^(X3)$^(X2) - $«(X3)«J>;(X2))] 

xi?(f(tf )42)(xi,X2,X3,6l,63)} (41) 



M'Te4 = 327rCFM|V2iVc Jo rfxicZx2dx3 /o°° 61^6163 

db^^BMM) X r;{[x3($^(x3) + $^(X3))$^(X2) 
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XEi''^'{t^J^)h'j\xi,X2,Xs,h,h) 

+ (l-X3)($;;(x3) + $^(ar3))$^(x2) 
X4'>(if (xi,X2,X3,6i,63)} (42) 



X^gg = -167rCFM|x/2iVc/o dx^dx^dx^ hdbi 

bsdh^BM, h) X r;{x3[$^(X3) - $»(X3)]$J(X2) 
X£;^«)'(4'^)4')(X1,X2,X3,61,63) 

+[(1 - X3) {^;{xs) - $«(X3)) $J(X2) - 2r; 

(1 - X3 + X2) ($;^(X3)$;(X2) - $^(X3)^>;(X2))] 

x4'^)'(if )/if (xi, X2, X3, 6i, 63)} (43) 

The evolution factors are given by 

E^^'{t)=as{t)af{t)Sm,=,, (44) 

with the Sudakov factor S = SBsS(j>Sp and the Wilson coefficients are given by 

a'l = C2/N,, 
< = 3/2C10/N,, 

a'e = 3/2Cs/Nc. (45) 

2.3 Numerical analysis 

The parameters used in our calculations are: the Fermi couphng constant Gp = 1.16639 x 10~^GeV~'^, 
the meson masses Mb = 5.28GeV, Mbs = 5.37GeF, Mp = 0.77GeV, = IMGeV, the decay 
constant fp = 0.205GeV, fj = 0.155Gey, = 0.237^6^", /J = Q.22QGeV , the central value of the 
CKM matrix elements 7 = 60°, \Vtd\ = 0.0075, \Vtb\ = 0.9992, = 0.0047, \Vus\ = 0.2196 and the 
meson lifetimer^ = l.Qbps,TBs — 1.461ps[?]. 

Using the above parameters, we get the branching ratios and helicity amplitudes of <pp 
decays (the helicity amplitudes are in Table 1) 

Br{B^ (j)p'^) = 4.1 X 10-^ 
Br(B° ^ (/)/9°) = 1.9 X 10-^ 
Br{B'^ <j)p°) = 3.09 X 10"^, 

Br{Bs° (t)p°) = 3.60 X 10-'^ (46) 

Compared with the averaged results of QCDF[?] 

BR{B- p^(t)) = 5.5 X 10-^ 
BR{B^ (/)p°) = 2.5 X 10-9. (47) 

or those from naive factorization [?] 

Bi?(5° ^ p^(j)) = 2.92 X 10-^ (48) 
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our predictions for B ^ (f)p are consistent with those of QCDF, the predictions for S° p^cf) are 
consistent with the result of naive factorization, because the nonfactorization effects in Bg —>■ p^(f) are 
httle. 

Presently, only the experimental upper limits are available at the 90% confidence level [?] , 

^ < 1.6 X lo-^ 

Br{B^ < 1.3 X 10-^ 

5^(5° ^ < 6.17 X 10"^. (49) 

Obviously, our results are consistent with the data. Our predictions will be tested by the oncoming 
measurements. 

For B ^ ^p decays, only penguin operators contribute, so there is no direct CP violation in the 
B ^ (f)p decays. 

For Bg —>■ (f)p^ decays, the CP asymmetry is time dependent 

Acp{t) ~ A^^p cos (Ami) + Ue+e' sin (Ami), (50) 

where Am is the mass difference of the two mass eigenstates of neutral Bg mesons. 
The direct CP violation parameter is defined 

- \M\^+\Mr ^ ' 

The direct CP violation parameter we can get in Bg (pp^ is 

^Cp{Bl^4>P^) =-8.0%. (52) 

From Table 1, we can find the longitudinal fraction in B^ cpp^ decays go beyond 70%, but the 
longitudinal fractions m. B ^ (f)p decays are very small, which is similar to B^ — p^p^ decay mode[?]. 
Id. B ^ (f)p decays, Oi,2iii 'Heff don't contribute via factorizable diagrams, the penguin operators 
contributing via factorizable diagrams are color suppressed, so that the nonfactorizable effects are 
the same order as the factorizable ones, which cause the B ^ cpp decays not to be factorization 
dominated, besides the above reason, the nonfactorizable longitudinal amplitude is opposite in sign 
to that of the factorizable part, therefore the longitudinal amplitude gets a large cancellation between 
the factorizable effects and the non factorizable parts such that \Aq\'^ is reduced much. 

According to Ref.[?], we can get the p,(f) vector meson transition form factors , 

B^ p, 1/(0) = 0.303, Ao (0) = 0.308, 

^i(O) = 0.233,^2(0) = 0.208, 
Bg 4>, V{0) = 0.430, Ao(0) = 0.363, 

Ai(0) = 0.304, yl2(0) = 0.276, 

which are consistent with the results with light cone sum rules [?] . 



3 Summary 

In this paper, we calculate the branching ratios and polarization fractions of 5 — ^ and Bg — > 
(j)p^ decays in perturbativc QCD approach, the predicted branching ratios are compared with the 
experimental data and results obtained with other approaches. CP parameters in Bg (pp^ are given 
in our paper, we compared with the experimental values, our results are consistent with the current 
experimental data. 
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